Formation of Cooper pairs between conduction and localized electrons in 

heavy-fermion superconductors 
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Cooper pairing between a conduction electron (c electron) and an / electron, referred to as the 
"c-f pairing," is examined to explain s-wave superconductivity in heavy-fermion systems. We first 
apply the Schrieffer- Wolff transformation to the periodic Anderson model assuming deep / level and 
strong Coulomb repulsion. The resulting effective Hamiltonian contains direct and spin-exchange 
interactions between c and / electrons, which are responsible for the formation of the c-f Cooper 
pairs. The mean-field analysis shows that the fully gapped c-f pairing phase with anisotropic s-wave 
symmetry appears in a large region of the phase diagram. We also find two different types of exotic 
c- f pairing phases, the Fulde-Ferrell and breached pairing phases. The formation of the c-f Cooper 
pairs is attributed to the fact that the strong Coulomb repulsion makes a quasiparticle / band near 
the center of the conduction band. 

PACS numbers: 74.70.Tx, 74.20.Mn, 74.25.Dw 



I. INTRODUCTION 

Various types of heavy-fermion superconductors dis- 
covered recently have attracted growing attention due 
to their unconventional features. Some materials with- 
out inversion symmetry have a superconducting phase in 
which the mixing of spin-singlet and spin-triplet states 
is expectedfi~— It has also been found that a multilayer 
material shows a strong-coupling superconducting state 
where the ratio of the superconducting gap to the tran- 
sition temperature, 2/S./k^T Cl is quite large compared 
to the conventional BCS valued Furthermore, possible 
signatures of the Fulde-Ferrell-Larkin-Ovchinnikov states 
have been observed in CeCoIns. 6 Although many differ- 
ent heavy-fermion superconductors have been found, the 
theoretical studies are still insufficient to deeply under- 
stand the individual superconducting properties. 

To reveal the mechanism of different types of supercon- 
ductivity, first of all, identifying the pairing symmetry is 
of crucial importance. In usual heavy-fermion supercon- 
ductors, the strong Coulomb repulsion between / elec- 
trons favors the nodal d-wave symmetry, which has been 
the subject of a number of theoretical studies, including 
the slave-boson approximation with l/A^-expansion^^ 
the random-phase approximation's^ the fluctuation- 
exchange approximation, 11 and the third-order pertur- 
bation approaches^ - — The experimental results also 
support the d-wave symmetry. The nuclear magnetic 
and quadrupole resonances (NMR and NQR) in typi- 
cal heavy-fermion compounds show a power-law temper- 
ature dependence of the spin-lattice relaxation rate and 
the lack of the Hebel-Slichter peak t 15 i 16 which indicate 
the existence of line nodes. Moreover, the phase dia- 
gram has the same feature as that of high-T c cuprates 
with ci-wave symmetry; superconductivity appears near 
the antiferromagnetic phased 

However, conventional s-wave superconductivity has 
also been found in some compounds. In NQR measure- 



ments on CeRu?i£ and CeCo2r^ the spin- lattice relax- 
ation rate exhibits an exponential decay at low tempera- 
tures and shows the Hebel-Slichter peak. Moreover, the 
recent photoemission spectroscopy (PES) experiment on 
CeRu^ has shown that the density of states (DOS) has 
a clear superconducting gap at the Fermi level. All these 
results were interpreted as evidence for the fully gapped 
pairing state with s-wave symmetry. Usually, this type of 
simple pairing symmetry can be understood within the 
framework of the conventional electron-phonon mecha- 
nism. However, it is unclear whether the electron-phonon 
attraction can be dominant since the Coulomb repulsion 
is rather strong in heavy-fermion systems. 

In this paper, we propose another possible way to un- 
derstand s-wave superconductivity in heavy-fermion sys- 
tems. The essence of our idea is to consider the Cooper 
pairing between a conduction electron (c electron) and a 
localized / electron, which we call the 11 c-f pairing." This 
type of Cooper pairing was previously examined in the 
study based on a slave-boson approach^ In this study, 
since the constraints on the enlarged Hilbert space are 
treated at the mean-field level, the effects of unphysical 
states are included in the solution. Using another the- 
oretical treatment, we present a detailed analysis of the 
c-f pairing state, including the derivation of the phase 
diagram, from a different point of view. By performing 
the Schrieffer- Wolff transformation to the periodic An- 
derson model, we first derive an effective Hamiltonian 
for deep / level and strong Coulomb repulsion. The re- 
sulting effective Hamiltonian includes direct and spin- 
exchange interactions between c and / electrons, which 
lead to the formation of the c-f Cooper pairs. We analyze 
the effective Hamiltonian within the mean-field approxi- 
mation and obtain the phase diagrams involving several 
types of c-f superconducting phases. Especially, we find 
the fully gapped state with anisotropic s-wave symmetry 
in a large region of the phase diagram. We also show 
that more exotic c-f pairing phases, the Fulde-Ferrell 
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(FF) and breached pairing (BP) phases, can appear in 
the other regions of the phase diagram. 

This paper is organized as follows. In Sec. UH we in- 
troduce the periodic Anderson model and derive an effec- 
tive Hamiltonian by using the Schrieffer- Wolff transfor- 
mation. We obtain the self-consistent equations for the 
order parameter of the c-f pairing superconductivity and 
some other quantities within the mean-field approxima- 
tion. In Sec. nn we show the results of our numerical 
calculations. We find three different types of c-f pairing 
phases in the ground-state phase diagram. At the end 
of the section, we discuss the reason for the formation of 
those c-f pairing states. Finally, Sec. IIVI is devoted to 
conclusions. 



II. MODEL AND CALCULATIONS 

We consider a typical heavy-fermion system composed 
of itinerant c electrons and nearly localized / electrons, 
which hybridize with each other. Usually, such a sys- 
tem is modeled by the periodic Anderson Hamiltonian 
-ffpAM — H + Hy> 

H ° = ~ l E 5>*W + H - c -) + E n L 

(ij) a ia 



H v = vJ2(fl^ + a-c), 



(2) 



where c\ a (fj a ) is the creation operator of a c electron 



(an / electron) with spin a at site i, n\ a = c\ a Ci a , and 
n la = ftrfiv Here, t is the hopping integral of c elec- 
trons, e/ is the position of the bare / level, p is the 
chemical potential, U is the on-site Coulomb repulsion in 
the / orbital, and V is the hybridization matrix element 
between c and / states. The sum (ij) runs over nearest- 
neighbor pairs of lattice sites. We consider the case of a 
square lattice in this study. 

In order to obtain an effective Hamiltonian describ- 
ing the c-f pairing superconductivity, we perform the 
Schrieffer- Wolff transformation 22 H = e s Hpamb~ s ', 
where S is chosen so as to eliminate all first-order terms in 
V. The generator S must satisfy the condition [S, H ] = 
—Hy, and is given by 



S = 



1 



E 



Ve 



-ik-Ri 



kv'ir 
y e — ik-R, 

£k-e/ 



U 



f t f 



(l - nCj cij ia H, 



(3) 



where a =f (J,) for a =\. (j), N is the total number 
of lattice sites, and e k = —2t(cosk x + cask y ). Here we 
set the lattice constant a = 1. When |e/| and €f + U 
are large compared to the effective kinetic energy of / 
electrons, which is roughly proportional to pV 2 (p is the 



c-electron DOS at the Fermi level), the system is approx- 
imated by keeping only the zeroth and second orders in 
V as HnHo + H 2 =H eS : 

H 2 = — [S, Hy] = -ffdir + ffcx + H c h + H p h , (4) 



where 



Hd " = ]j E (^ k ' k - J J k'k(nf t + nf ; ) 

k'kicr 



h, 



xe - i(k '- k ).R ic t^ CkCT) (5) 



k'ki 

+ S r c l't c ^ + ^ 2 ( c k't ck t~ c k'i c ki)) , (6) 
ff ch = -^E (^ kk - 1 Jkk ^ + n ^)) 

Uija 

Ke -ik.(R i - %)/ j CT/ . (T) (7) 

H P h = iA? E Jk ' k ' 

x( e - l(k ' +k)R '4, ff cL/ JCT / i5 +H.c.). (8) 

Here, Sj = ^J2a>afL' T <?'<?f™ is the s P in operator of / 
electrons, and the coupling energies Jk'k and Wk'k are 
defined as 



V 

Jk'k = ~rj~ (ik' + Lu) , 

WW = ^ {M v + M k ) , 



(9) 
(10) 



with L k = f/ 2 (e k — e/) 1 (ek — e/ — U) 1 and Af k 
I ' (e k — e/)~\ As shown in Eqs. (TJJ-flH]), the second term 
Hi consists of four different interaction terms: the direct 
interaction ffdir, the spin-exchange interaction _ff C x, the 
/-electron correlated hopping iJ c h, and the pair hopping 

We apply the mean-field approximation to many-body 
terms in H c r . We first introduce the following order pa- 
rameter characterizing c-f superconducting phases: 



Ak= — E ^'k^k', 



with 



(fl 



k'+qt -k'4- 



/-k'+ql C k't) 



(11) 



(12) 



By decoupling Hdir and -Hex, one can extract the c-f su- 
perconducting order parameters, which means that these 
terms play a crucial role for the formation of the c-f 
Cooper pairs. The effective mass of c electrons is much 
smaller than that of / electrons. Recently, this type of 
Cooper pairing with unequal masses has been intensively 
studied in the field of ultracold Fermi gasesi 2 ^— Liu and 
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Wilczek have discussed this issue by assuming an attrac- 
tive interaction between fermions with different masses^ 
They found that the mean-field phase diagram contains 
two different types of unconventional superfluid phases as 
well as the usual fully gapped s-wave superfluid phase. 
One of them is the FF phase with a finite center-of-mass 
momentum of the Cooper pairs^ and the other is the 
BP phase^Sr— which was called the interior gap super- 
fluid phase in the original paper. In the BP phase, the 
Cooper pairs have zero center-of-mass momentum, while 
the Bogoliubov band has no gap unlike the case of the 
fully gapped s-wave state. The name "breached pair- 
ing" comes from the fact that the superfluid component 
is "breached" by the normal fluid component.— To take 
into account the possibility of the FF state, we assume 
a finite center-of-mass momentum q of the Cooper pairs 
in Eq. (ff"2]) . The order parameter A k can be chosen to 
be real without loss of generality. 

In addition to the superconducting order parameter, 
we also include all possible Hartree-type mean fields, 
which are defined as 



k 
k 

c =jjJ2Jkk(ci a c ka ), a =1,1, 

k 



(13) 
(14) 
(15) 

(16) 



where n=n c + n/ is the total density of the system. De- 
coupling the Hubbard term in Hq and all the terms in Hi , 
we obtain the following mean- field Hamiltonian Hmf'- 

H M F = E 6<C k(T C k(T + 4/L/k<r 

+ 2 E A k ( c kt/-k+ q i - /k +q /-ki + H.c.) 

i£A k fl k + iwn/6 (17) 



1 
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f 



where 



; k+qJ 



A 
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(19) 
(20) 



and f(E) = l/(e^ E +l), with^ = 1/Tis the Fermi distri- 
bution function. The upper (w k ) and lower (w k ) Bogoli- 
ubov bands are related to £ k by cj k = max |, |£ k |) 
= min(|£+|,|£ k |). The conditions = 0, 



and 



i an dn 

iV 3/i' dcp 



= 0, and — give the self-consistent 

1 on f ° 



equations for A k , n, n/, and <f>, respectively: 



A k = 



1 V- j A /(^)-/(^k") 



C k'+q 
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E 
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(21) 
(22) 



L E 
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(23) 



2N 



(Jk'+q + Jk0(/(^)+/(^k")) 



(Jk'+q — ^k')(^k' 



3 k'+q 



Ekf 

x(f(£+)-f(£ k ,)) 



(24) 



with J k = J kk . We can see from Eq. © that the or- 
der parameter A k can be separated into constant and 
k-dependent parts as A k = Ao + AiL k . Substituting 
this expression, we derive the equations for Ao and Ai 
instead of Eq. (j2"Tj) : 



with £ k = e k -/x + W kk -n/J kk /4, e[ = e f -fi + Un f /2- 
(f)/2 - Wkk + n/Jkk/4, and <j> = 4> c - 4> f . Note that in 
Eq. (|17|) . the effective one-body energy of / electrons, 
e£, depends on the wave vector k. This means that the 
correlations between c and / states yield a finite band- 
width for / electrons. From Eq. (|17|) . we easily find the 
corresponding thermodynamic potential 



^=E 
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-A k B k 
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(18) 



k,a=± 



Ao = 77^E ik '( Ao + A i Lk 



UN 
V 



A i-77^E( A o + A i L k' 



f(S+)-f(Sy) 



(25) 



(26) 



For a given total density n, the values of Ao, Ai, /i, rz/, 
and 4> are obtained by numerically solving Eqs. (|22p - (|2"6")l 
in a self-consistent way. At the same time, we also need 
to minimize the free energy F = fi + finN with respect to 
q. In the present study, the center-of-mass momentum 
q is assumed as q = q/iJ2), and the Coulomb 

repulsion U is set to be U/t — 12. 
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FIG. 1. The phase diagrams in the (|e/|/i, V/t) plane at 
n = 2.2 and T/i = 0.005. The solid curves indicate the 
second-order phase transitions between the fully gapped s 
wave, FF, BP, and normal phases. The dashed curves indi- 
cate the position where the gap in the lower Bogoliubov band 
of the FF phase vanishes. The dotted vertical lines represent 
the lines of \tf\/t = 5.5 and 5.8, respectively. 
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FIG. 2. The V/t dependences of \A \/t, \Ai\/t, q, and A/t 
at n = 2.2, T/t = 0.005, and \e f \/t = 5.8. The solid vertical 
lines indicate the second-order phase transitions. The dashed 
vertical line represents the position where the gap A of the 
FF phase vanishes. 



III. RESULTS 

Before presenting the results of the calculations, we 
briefly comment on the difference between our study and 
the previous work by Hanzawa and Yosida. 21 Hanzawa 
and Yosida discussed the c-f pairing state on the basis 
of the periodic Anderson model in the limit of strong 
Coulomb repulsion, where the doubly occupied states in 
the / orbital are excluded from the Hilbert space. They 
derived the gap equation and estimated the order of the 
transition temperature for the c-f pairing superconduc- 
tivity. In our present study, we further take into account 
the effect of Hartree-type mean fields and the possibility 
of the pairing state with finite center-of-mass momen- 



tum. Since the Coulomb repulsion is large but finite in 
our analysis, the influence of doubly occupied states is in- 
cluded in the results. As a consequence, we find several 
unconventional c-f pairing phases in addition to the sim- 
ple c-f pairing phase discussed by Hanzawa and Yosida. 

Figure [T] shows the |e/|-V phase diagram at n = 2.2, 
which includes four different phases: the fully gapped s 
wave, FF, BP, and normal phases. Let us discuss the 
phase transitions between these phases along the line of 
\tf\/t — 5.8, which is depicted by the dotted vertical line 
in Fig. [T] We show the V/t dependencies of A , Ai, 
and q at \cf\/t = 5.8 in Fig. [5J We also show the ac- 
tual gap in the lower Bogoliubov band, A=mink (w^"). 
Since the sign of Ao is always opposite to that of Ai 
in the parameter range of Fig. [2j we plotted the ab- 
solute values |A | and |Ai| in the figure. For large V, 
the fully gapped s-wave phase is preferred. The order 
parameter Ak has anisotropic s-wave symmetry and the 
corresponding lower Bogoliubov band u;^ shows a finite 
gap, as shown in Figs. [3][a) and[3][b). As V decreases, 
the FF state appears as the ground state. Due to the 
existence of finite q, the lower Bogoliubov band has 
an asymmetry with respect to the center of the Brillouin 
zone, as seen in Fig. SJa). The band is gapless, namely, 
A = 0, in most of the FF region. Only in a narrow re- 
gion (1.508<V/t<1.524) of Fig. H]do we have the fully 
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FIG. 3. (a) The absolute value of the superconducting order 
parameter |Ak| and (b) the lower Bogoliubov band ui^ in the 
fully gapped s-wave state at n = 2.2, T/t = 0.005, |e/|/t = 
5.8, and V/t = 1.8. 
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FIG. 4. The lower Bogoliubov bands at n = 2.2, T/t - 
0.005, and \e f \/t = 5.8 for (a) V/t = 1.3 and (b) V/t = 1, 
respectively. 



gapped FF state. As V is decreased further, the tran- 
sition to the BP phase occurs at V/twl.025, where the 
center-of-mass momentum q vanishes. As shown in Fig. 
IDJb), the lower Bogoliubov band uj^ of this phase touches 
the zero-energy line, although it is symmetric about the 
center of the Brillouin zone. For even smaller V, we have 
only the trivial solution Ak = 0, which is natural since 
the effective attraction | Jkk' | between c and / electrons 
[Eq. (|10p ] becomes smaller as V is decreased. 

It is worthy to note that the |e/|-F phase diagrams for 
n = 2 + S and n = 2 — S are symmetric with each other 
about |e/| = U/2 = 6t. For example, we can obtain the 
phase diagram for n = 1.8 by the left-right inversion of 
Fig. [T] This symmetric property comes from the fact 
that the periodic Anderson model has particle-hole sym- 
metry at e/ = — U/2 in the case of bipartite lattices. 

Next, let us examine the effect of temperature on the o 
/ pairing phases. Figures [5ja) and [5jb) show the phase 
diagrams in the (T/t, V/t) plane for \ef\/t = 5.5 and 
\cf\/t = 5.8, marked by the dotted vertical lines in Fig. [T] 
We can see that the fully gapped s-wave pairing state is 
more robust against temperature than the nodal pairing 
states. Especially, the FF phase completely disappears 
as the temperature is increased. Such a sensitive tem- 
perature dependence of the FF phase has been obtained 
in previous studies! 32 ' 33 The region of the BP phase also 
gets smaller with increasing temperature, but it still sur- 
vives after the disappearance of the FF phase. For high 
temperatures, the fully gapped s-wave phase occupies a 
large region of the phase diagram. 

Since the bare / level e/ is rather deep below the Fermi 
energy, we should clarify the reason why a c electron 
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order phase transitions between the fully gapped s wave, FF, 
BP, and normal phases. The dashed curves indicate the po- 
sition where the gap A of the FF phase vanishes. 



near the Fermi level can form a pair with an / electron. 
We start the discussion from the original periodic An- 
derson model given by Eqs. ([T]) and ([2]). Because of 
the hybridization V between c and / states, the effective 
/ level has a finite dispersion, namely, a finite effective 
mass, even if the bare / band is completely flat. Fur- 
thermore, the strong Coulomb repulsion U splits the ef- 
fective / level into the upper and lower Hubbard bands 
and forms a quasiparticle / band in between them, as in 
the case of the standard Hubbard model^ The forma- 
tion of the quasiparticle / band has been shown by the 
previous studies using the dynamical mean-field theory 
(DMFT)j24"^S Especially near the half- filling, the quasi- 
particle / band is generated in the vicinity of the Fermi 
level, i.e., near the center of the conduction band. This 
allows us to propose that a conduction electron forms a 
pair with an electron in the quasiparticle / band and it 
causes superconductivity. 

In the present work, the effect of the Coulomb repul- 
sion U is treated within the mean-field approximation, in 
which the splitting of the effective / level is not described. 
However, the quasiparticle / band is approximately ex- 
pressed by the Hartree shift as e/ = e/ + Unf/2. Fig- 
ureslUJa) and[6jb) show an example of the partial DOSs 
(PDOSs) of the c and / bands in the fully gapped s-wave 
state. The PDOS of / electrons has a large weight near 
the Fermi level and the superconducting gap opens in 
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the fully gapped s-wave state at n — 2.2, T/t — 0.005, |e/|/t = 
5.8, and V/t = 1.8. The dashed vertical lines represent the 
Fermi level. 



both the PDOSs, which support our scenario proposed 
above. It should be noted, however, that our mean-field 
treatment may overestimate the PDOS of the effective / 
level near the Fermi energy. 

It is known that at half-filling the periodic Anderson 
model has an insulating ground stated - — which exhibits 
antiferromagnetic order when the Coulomb repulsion is 
larger than the critical value [7 c 42r— In the case of fi- 
nite doping, the self-consistent second-order perturbation 
approach by Mutou^ showed that this model favors a 



metallic ground state, in which the quasiparticle / band 
is located around the Fermi level. However, this study 
did not take into account the c-f pairing state. We ex- 
pect that the c-f pairing state can appear in a doped 
region of the periodic Anderson model. In order to dis- 
cuss the doping-induced phase transition from insulator 
to c-f pairing state, it is required to perform a further 
analysis which can treat the insulating states, e.g., the 
use of the DMFT, although it is beyond the scope of the 
present study. 



IV. CONCLUSION 

We have studied the Cooper pairing between a conduc- 
tion electron (c electron) and an / electron, called the "c- 
/ pairing,"— to understand s-wave superconductivity in 
heavy-fermion systems. Considering a system with deep 
/ level and strong Coulomb repulsion, we first derived 
an effective Hamiltonian by performing the Schrieffer- 
Wolff transformation to the periodic Anderson model. 
Within the mean-field analysis of the effective Hamilto- 
nian, we obtained the ground-state phase diagrams in- 
cluding three different types of c-f pairing phases: the 
fully gapped, FF, and BP phases. Especially, we found 
that the fully gapped c-f pairing phase with anisotropic 
s-wave symmetry occupies a large region of the phase di- 
agram. Moreover, we demonstrated that the fully gapped 
c-f pairing state is more robust against temperature than 
the FF and BP phases. Our results may be relevant to 
the recent experiment which observed an anisotropic s- 
wave superconducting gap in CeRu2^ 



ACKNOWLEDGMENTS 

We thank Grant-in- Aid from JSPS (K. M.) and KAK- 
ENHI (23840054) from JSPS (D. Y.). 



* masuda@kh.phys.waseda.ac.jp 

1 E. Bauer, G. Hilscher, H. Michor, Ch. Paul, E. W. Scheidt, 
A. Gribanov, Y. Seropegin, H. Noel, M. Sigrist, and P. 
Rogl, Phys. Rev. Lett. 92, 027003 (2004). 

2 N. Kimura, K. Ito, K. Saitoh, Y. Umeda, H. Aoki, and T. 
Terashima, Phys. Rev. Lett. 95, 247004 (2005). 

3 V. M. Edelstein, Sov. Phys. JETP 68, 1244 (1989). 

4 L. P. Gor'kov and E. I. Rashba, Phys. Rev. Lett. 87, 
037004 (2001). 

5 Y. Mizukami, H. Shishido, T. Shibauchi, M. Shimozawa, 
S. Yasumoto, D. Watanabe, M. Yamashita, H. Ikeda, T. 
Terashima, H. Kontani, and Y. Matsuda, Nature Phys. 7, 
849 (2011). 

6 Y. Matsuda and H. Shimahara, J. Phys. Soc. Jpn. 76, 
051005 (2007). 

7 M. Lavagna, A. J. Millis, and P. A. Lee, Phys. Rev. Lett. 
58, 266 (1987). 



A. Houghton, N. Read, and H. Won, Phys. Rev. B 37, 
3782 (1988). 

D. J. Scalapino, E. Loh, and J. E. Hirsch, Phys. Rev. B 
34, 8190 (1986). 

K. Miyake, S. Schmitt-Rink, and C. M. Varma, Phys. Rev. 
B 34, 6554 (1986). 

R. Arita, K. Kuroki, and H. Aoki, Phys. Rev. B 60, 14585 
(1999). 

H. Ikeda, J. Phys. Soc. Jpn. 71, 1126 (2002). 

Y. Nishikawa, H. Ikeda, and K. Yamada, J. Phys. Soc. Jpn. 

71, 1140 (2002). 

H. Fukazawa and K. Yamada, J. Phys. Soc. Jpn. 72, 2449 
(2003). 

T. Mito, S. Kawasaki, G.-q. Zheng, Y. Kawasaki, K. Ishida, 
Y. Kitaoka, D. Aoki, Y. Haga, and Y. Onuki, Phys. Rev. 
B 63, 220507(R) (2001). 

Y. Kohori, Y. Yamato, Y. Iwamoto, T. Kohara, E. D. 
Bauer, M. B. Maple, and J. L. Sarrao, Phys. Rev. B 64, 



7 



17 



134526 (2001). 

C. Pfleiderer, Rev. Mod. Phys. 81, 1551 (2009). 



18 K. Matsuda, Y. Kohori, and T. Kohara, J. Phys. Soc. Jpn. 
64, 2750 (1995). 

19 K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, H. Sug- 
awara, Y. Aoki, and H. Sato, Physica B 237, 304 (1997). 

20 T. Kiss, F. Kanetaka, T. Yokoya, T. Shimojima, K. Kanai, 
S. Shin, Y. Onuki, T. Togashi, C. Zhang, C. T. Chen, and 
S. Watanabe, Phys. Rev. Lett. 94, 057001 (2005). 

21 K. Hanzawa and K. Yosida, J. Phys. Soc. Jpn. 56, 3440 
(1987). 

22 P. Sinjukow and W. Nolting, Phys. Rev. B 65, 212303 

(2002) . 

23 W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 

(2003) . 

24 M. Iskin and C. A. R. Sa de Melo, Phys. Rev. Lett. 97, 
100404 (2006). 

25 G.-D. Lin, W. Yi, and L.-M. Duan, Phys. Rev. A 74, 
031604(R) (2006). 

26 S.-T. Wu, C.-H. Pao, and S.-K. Yip, Phys. Rev. B 74, 
224504 (2006). 

27 M. M. Parish, F. M. Marchetti, A. Lamacraft, and B. D. 
Simons, Phys. Rev. Lett. 98, 160402 (2007). 

28 P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 

29 E. Gubankova, W. V. Liu, and F. Wilczek, Phys. Rev. 
Lett. 91, 032001 (2003). 

30 Michael McNeil Forbes, E. Gubankova, W. V. Liu, and F. 
Wilczek, Phys. Rev. Lett. 94, 017001 (2005). 

31 L. He, M. Jin, and P. Zhuang, Phys. Rev. B 73, 214527 
(2006). 



32 T. K. Koponen, T. Paananen, J. -P. Martikainen, and P. 
Torma, Phys. Rev. Lett. 99, 120403 (2007). 

33 T. K. Koponen, T. Paananen, J.-P. Martikainen, M. R. 
Bakhtiari, and P. Torma, New J. Phys. 10, 045014 (2008). 

34 A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

35 Y. Shimizu, O. Sakai, and A. C. Hewson, J. Phys. Soc. 
Jpn. 69, 1777 (2000). 

36 L. de' Medici, A. Georges, G. Kotliar, and S. Biermann, 
Phys. Rev. Lett. 95, 066402 (2005). 

37 G. Sordi, A. Amaricci, and M. J. Rozenberg, Phys. Rev. 
Lett. 99, 196403 (2007). 

38 A. Amaricci, G. Sordi, and M. J. Rozenberg, Phys. Rev. 
Lett. 101, 146403 (2008). 

39 M. Jarrell, H. Akhlaghpour, and Th. Pruschke, Phys. Rev. 
Lett. 70, 1670 (1993). 

40 T. Mutou and D. S. Hirashima, J. Phys. Soc. Jpn. 63, 4475 
(1994). 

41 T. Pruschke, R. Bulla, and M. Jarrell, Phys. Rev. B 61, 
12799 (2000). 

42 M. J. Rozenberg, Phys. Rev. B 52, 7369 (1995). 

43 M. Vekic, J. W. Cannon, D. J. Scalapino, R. T. Scalettar, 
and R. L. Sugar, Phys. Rev. Lett. 74, 2367 (1995). 

44 S. Horiuchi, S. Kudo, T. Shirakawa, and Y. Ohta, Phys. 
Rev. B 78, 155128 (2008). 

45 T. Mutou, Phys. Rev. B 62, 15 589 (2000). 



